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Abstract 



The lifetimes of metastable states in kinetic Ising ferromagnets are studied 
f^ ■ by droplet theory and Monte Carlo simulation, in order to determine their 

^ . dependences on applied field and system size. For a wide range of fields, the 

dominant field dependence is universal for local dynamics and has the form 
C^ I of an exponential in the inverse field, modified by universal and nonuniversal 

Vh ' multiplicative power-law prefactors. Quantitative droplet-theory predictions 

^-^ ■ for these dependences are numerically verified, and small deviations from the 

G . predictions are shown to depend nonuniversally on the details of the dynam- 

Q i ics. We identify four distinct field intervals in which the field dependence 

and statistical properties of the lifetimes are markedly different. The field 
L^ marking the crossover between the weak-field regime, in which the decay is 

j_j ■ dominated by a single droplet, and the intermediate- field regime, in which it 

is dominated by a finite density of droplets, vanishes logarithmically with sys- 
tem size. As a consequence the slow decay characteristic of the former regime 
may be observable in systems that are macroscopic as far as their equilibrium 
properties are concerned. 
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I. INTRODUCTION 



Although it is observed in nature in contexts as dissimilar as supercooled water or water 
vapor |jl|-0] and the electroweak |^ and QCD confinement/deconfinement [Q phase transi- 
tions, met ast ability is very difficult to characterize in a microscopically precise fashion . In 
certain systems with weak long-range interactions infinitely long-lived metastable states can 
exist in the thermodynamic limit 0, and metastability in such systems has been studied with 
several techniques [Q, including field theoretical, Monte Carlo (MC), and transfer-matrix 
methods. However, in systems with short-range interactions metastable states eventually 
decay, even though their lifetimes may be many orders of magnitude larger than other char- 
acteristic timescales of the system, and may even become comparable to the age of the 
universe [Q. In order to gain a deeper understanding of metastability in short-range sys- 
tems, it is necessary to consider in detail the lifetimes of metastable states and how they are 
determined by the physical mechanisms involved in the decay [Q. 

As a prototype for the metastable dynamics of short-range systems, the decay of the 
magnetization in impurity-free kinetic Ising ferromagnets in unfavorable applied fields has 
been studied by MC methods in both two |TO|-p!^] and three |l9|-pl|] dimensions. The results 
of several of these studies were analysed in terms of droplet theory |]22|-p8|, establishing gen- 
eral agreement between theory and simulations. In the present work we further investigate 
the extent of that agreement by analytical droplet-theory calculations and MC simulations. 
We emphasize the dependences on applied field and system size of the metastable hfetimes 
and their relations to particular decay mechanisms. Our analysis extends a recent study by 
two of us [pQJ^ , and Ref. |2^ is henceforth referred to as (I). 



The main results in (I) can be summarized as follows. In agreement with previous |n,20 



and more recent pl| simulations it was found that the statistical properties and system-size 
dependence of the metastable lifetimes in MC simulations were markedly different in two 
separate field regions. The point separating these regions was called "the dynamic spinodal 
point" (DSP), and the corresponding field, which depends on temperature and system size, 
was denoted -ffosp- For \H\ > i^DSP, the mean lifetime was observed to be independent 
of system size and much greater than its standard deviation [0,^,^. This field region 
was therefore called "the deterministic region". For weaker fields, the mean lifetime was 



inversely proportional to the system volume |I^, p!3| , P(]| , pT| , P^ -p5| , and the standard deviation 
was approximately equal to the mean |[T^j20|j21[] . For these reasons this regime was termed 
"the stochastic region" . The position of the DSP, which separates these regions, depends on 
the temperature and system size. The main goal of the present work is to understand and 
quantify these dependences. 

In this paper we present a quantitative analysis of metastable lifetimes in terms of a 
droplet model of homogeneous nucleation in c?- dimensional systems p2|-p^ , obtaining ana- 



lytic results for the logarithmic derivative of the lifetime with respect to 1/\H\'^~^, as well 
as the explicit size and temperature dependence of -^dsp- Both the deterministic and the 
stochastic region are found to be divided into two subregions, and expressions for the cor- 
responding crossover fields are derived. These analytic results are shown to be in excellent 
agreement with MC simulations for d=2 at a temperature of 0.8Tc. 

In order of increasingly strong unfavorable field the distinct regions that we identify are: 
the "coexistence region", characterized by subcritical fiuctuations on the scale of the system 



volume; the "single-droplet region", characterized by decay via a single critical droplet; 
the "multi-droplet region" , characterized by decay via a finite density of droplets; and the 
"strong-field region", in which the droplet picture is inappropriate. The two former regions 
comprise the stochastic region, and the two latter ones the deterministic region. The DSP 
thus marks the crossover between the single-droplet and multi-droplet regions. The crossover 
between the coexistence region and the single-droplet region was called "the thermodynamic 
spinodal point" (THSP) in (I), and the crossover between the strong-field and multi-droplet 
regions was called "the mean-field spinodal point" (MFSP). These different regions and the 
crossover points that separate them are summarized in Fig. |T], which shows a schematic 
sketch of the field dependence of the average metastable lifetime. 

Whereas the dominant field dependence of the lifetime in the single-droplet and multi- 
droplet regions is exponential in the inverse field, our MC data are sufficiently precise to also 
allow evaluation of multiplicative power-law prefactors. These are predicted by droplet the- 
ory P^|-P7[] but are difficult to detect numerically and, to our knowledge, have not previously 
been measured in MC studies of lifetimes or nucleation rates. 

The rest of the paper is organized as follows. In Sec. ^ the kinetic Ising model is 



defined, and the numerical methods used in this work are discussed. In Sec. |T| the droplet- 
theoretical predictions are developed, for both infinitely large and finite systems. In Sec. ^ 
the numerical results are presented and compared with the theory, and Sec. [V| contains a 
brief discussion. 



II. MODEL AND NUMERICAL METHODS 

The model is defined by the reduced Hamiltonian, 

nik^T = -kY. s^^j -hJ2s^^ (1) 

where Si = ±1 is the spin at site i, T is the temperature, ks is Boltzmann's constant, 
K = J/k-BT (>0) and H = h/k-QT are the reduced coupling constant and field, respectively. 
The sums Z](jj) and Y.i run over all nearest-neighbor pairs and over all A^ sites on a d- 
dimensional hypercubic lattice, respectively. 

The dynamics is given by the Metropolis single-spin-fiip MC algorithm. The transition 



probability for a fiip of the spin at site a from Sa to —Sa is thus defined as [^ 

W{sa -^ -Sa) = min[l,exp(-A£'Q,)] , (2a) 

where 



AEa = 2Sa 



NN 



(2b) 



is the energy change during the fiip, divided by k^T. 
We study the relaxation of the magnetization. 



N 

m{t)=N-'Y.s,it), (3) 



starting from an initial state magnetized opposite to the applied field. This approach was 
also used, e.g. in Ref. [^. The magnetization m{t) is directly obtained as an average 
over the droplet size distribution |r8| , pT| , p2[ , and it is closely related to the nonequilibrium 



relaxation functions introduced by Binder ||Tl[. The volume fraction of stable phase at time 



t is (j)s{t) = {mnis — 'm{t))/{m^s — 'ms), where rris and m^s are the bulk equilibrium and 
metastable magnetizations, respectively. 

The simulations were performed on square Lx L systems, with L=64, 128, 256, 400, 
and 720 and periodic boundary conditions. The systems were initialized with ??t,(0)=+1 and 
allowed to develop in a constant field H<0 at T = 0.8Tc (i^ = 0.550 858...). This temperature 
is high enough to obtain reasonable MC acceptance rates, but it is sufficiently low to avoid 
complications such as critical slowing-down and possible finite-size scaling effects due to 
small droplets [^]. For later reference we note that the equilibrium magnetization at 0.8Tc 



is ms = 0.954 410... |3^, and the equilibrium surface tension along a primitive lattice vector 



is (70 = 0.745 915... M 



Average first-passage times to m=0.9, 0.8, . . . , 0.0 were recorded, together with their 
empirical standard deviations, and the mean lifetime of the metastable state was estimated 
as the average first-passage time to m=0.7. We used this particular cutoff magnetization 
for the following reasons. In order for the droplet-theoretical considerations developed in 
Sec. |T| to remain valid during the whole time evolution until the cutoff, the cutoff value 
of 0s should be below the percolation limit. However, the cutoff magnetization should not 
be so close to mms that subcritical fiuctuations in the metastable phase ( "recrossing events" 
||T6| ) might be mistaken for decay events, thus leading to significant underestimation of the 
mean lifetime. The cutoff magnetization m=0.7, which, assuming m^s = -f^s^ corresponds 
to 0s ~ 0.13, was chosen as a reasonable compromise. However, except for very strong and 
very weak fields, the dependence of the observed lifetimes on the cutoff is weak, and results 
almost identical to those reported here were obtained with cutoff at m=0.0, or 0s ~ 0.5. 

A more accurate, but also more computationally intensive, method to estimate the life- 
times would be to use the recrossing events by choosing as a field-dependent cutoff a value 



of m which lies on the side of the maximum of the recrossing-event distribution ||T6[ oppo- 
site from that of the initial magnetization, and for which the recrossing-event probability is 
below a given threshold. We do not implement this method here, but plan to do so in future 
studies. 

For fields at which the lifetime is large, the magnetization m{t) was measured after each 
MC step per site (MCS), whereas for shorter lifetimes it was measured after each successfuU 
spin fiip. The latter method was used whenever necessary in order to ensure that the 
uncertainty due to the discreteness of the time variable was at least one order of magnitude 
smaller than the standard deviations in the average first-passage times. 



The numerical data were obtained with a special-purpose m-TIS2 computer P^J36|| at Ky- 
oto University and on heterogeneous clusters of IBM RS/6000, DEC-station, and DEC 3000 
Model 400 Alpha workstations at Florida State University. 

The m-TIS2 architecture only allows the MC updates to be performed sequentially, and 
the available memory limited the maximum system size to L=720. To assess the effects of the 
former restriction, simulations with updates at randomly chosen sites were performed on the 
workstation clusters for L=128 and 720. The total computer time spent was approximately 
400 m-TIS2 hours and 1600 hours of workstation time for the sequential-update study, and 



an additional 8000 workstation hours for the random-update investigation. 

III. DROPLET THEORY 

In this section we derive the droplet-theoretical relations necessary to interpret the sim- 
ulated lifetime data. 

The system is characterized by six length scales: the lattice constant, which we take as 
unity; the single-phase correlation lengths in the stable and metastable phases, ^s and .^ms, 
respectively; the critical droplet radius Re] the size Rq to which one droplet can grow before 
it is likely to meet another, which we shall refer to as the mean droplet separation; and 
the system size L. Here we only consider cases in which ^s is smaller than the other length 
scales, i.e. well below the critical temperature. Thus we are left to consider the interplay 
between four lengths: L, Rq, Rc, and ^ms- 

A. Infinite systems 

First we obtain i?o and -Rc for d-dimensional systems in the limit L — >■ cxd. Although 
the droplets are almost circular at our simulation temperature, we give a general argument 
that remains valid for low temperatures, at which nonspherical droplets appear because of 
the anisotropy of the surface tension [p^ , p^^ -|411] . The free energy of a (regular, but not 
necessarily spherical) (i-dimensional droplet of radius R (defined as half the extent of the 
droplet along a primitive lattice vector) and volume V{R)=QdR'^ is 

F{R) = \/('^-i)/'^S - VA 

= Q.^t^^/'^R'^-^t - QdR'^A , (4) 

where A is the difference in bulk free-energy density between the metastable and stable 
states. The quantity E is a temperature- and, in principle, field-dependent proportional- 
ity factor which relates the surface contribution to F{R) with the droplet volume |]3^ , ^ . 
Applying standard droplet-theory arguments |27] to F{R) one finds the critical radius, the 



free-energy cost of a critical droplet, and the nucleation rate per unit time and volume. The 
critical radius is 

p (^ TJ\ (^~^) ^ {d-l)ao {d-l)(ro 



A dn'J'' A 2mMT\H\ 



'd 



.1/(1 



where uo = T,/{dnJ ) is the surface tension along a primitive lattice vector p4|. The 



approximation Afii2mskBT\H\, where rUsiT) is the spontaneous equilibrium magnetization. 



is expected to be valid not only at low temperatures, but even near the critical point |^ 
The free-energy cost of a critical droplet is 



Fc{T,H) 



'd-lY'^ (t' 



A ; \d) 
d-i V-^ (t' 



2msk-BT\H\ 



(6) 



and the nucleation rate per unit time and volume, I{T,H), is determined by Fc through 



I{T,H) = A{T)\H\'+'e ~^^^ ^ A{T)\H\''~^''e m"-' 



with 



'd-iV-'f s V 



2ms J V dhsT 



(7a) 



(7b) 



where the approximation introduced in the second part of Eq. ( [7a| ) is the same as in Eqs. (^ 
and (^. The function A{T) is expected to be nonuniversal, 6 is a universal exponent related 
to Goldstone-mode "wobbles" on the droplet surfaces p3|,^, and c gives the H dependence 
of the "kinetic prefactor" p^j25[] , which is the only part of /(T, H) that may depend explicitly 
on the specific dynamics. A field-theoretical calculation gives [B^ 



{3-d)d/2 

-7/3 



for l<d<5, dj^3 
for d=3 



(8) 



Strictly speaking, a multiplicative correction term of form [1 + 0{H^)] should be included 
in the exponent in Eq. (|7aD [p^ , p^ - |41|j43|j4^ , but its effect on the metastable lifetimes is 
small for the relatively weak fields on which we focus our attention in this work, and it has 
therefore been suppressed for simplicity. 

For d=2 there is substantial numerical evidence that 6=1, as predicted by Eq. (§). This 
is obtained from calculations that do not involve the dynamics, such as analyses of series 
expansions pQJ^ , ^ and transfer-matrix calculations 
MC work 



3,P 



_^ ^,^. These studies, as well as 

also indicate that the free-energy cost of the critical droplet is given 
by Eq. (^) with the zero-field equilibrium value for E. We therefore adopt the notations 
S(T) and S(T) to emphasize the lack of field dependence in these quantities. They can 
be obtained with arbitrary numerical precision by combining a Wulff construction with 
the exact, anisotropic zero-field surface tension |5^. This general result, that the surface 
free energy of compact critical droplets is determined by the zero-field equilibrium surface 
tension, is also supported by MC studies of nucleation rates in three dimensions |T^-^. For 
dynamics that can be described by a Fokker-Planck equation it is expected that the kinetic 
prefactor is proportional to R~^ p^-p^, which by Eq. (|^) would yield c=2. However, we 
are not aware that independent, numerical verifications of this result have been performed 



previously. As discussed in Sec. ^^ we find that the measured value of c depends on the MC 
update scheme. 

Assuming that the growing droplets are not substantially deformed [0 , the radial growth 
velocity is obtained in an Allen-Cahn approximation p7|j^J^ as 



v± 



(rf-i)r (i?-^-i?-i) ^=^ (d-i)ri?-i = vo , 



(9) 



where F depends on the details of the kinetics. The mean droplet separation Rq and its 
associated time scale to are obtained from v± and the nucleation rate / by requiring that 
f j_to and RqIqI 



Rq - 
Eqs. 



1 |49|. If Rq ^ Re, so that v± ~ Vq, then by using vq oc \H\ from 



and 



one obtains 



b+c+d 



to{T,H) = (vp)-— = B{T)\H\-^s^ 6"+^ 



1 am 



(10a) 



and 



Ro{T,H) =voto = C{T)\H\-^^e^+'T^ 



(10b) 



The coefficients B and C have the same nonuniversal status as A in Eq. (|7a]) . Although 
i?c~ \H\~^ ^oo as \H\ — >-0, Rc/Rq^O in the same hmit. 

If one assumes that the positions of the (possibly overlapping) growing droplets are 
uncorrelated, the volume fraction of stable phase at time t becomes 



,(t) = 1 - exp 



n. 



d+1 \toJ 



d+r 



(11) 



which is known as Avrami's law P9|-|52|. From this it is seen that the average time taken to 
reach a specific 0s is 



(t(0s)) = to(r,ff) 



d+l 



n. 



ln(l- 



i/(d+i) 



(12) 



The "ideal-gas" approximation leading to Avrami's law is expected to hold when the total 
volume fraction of droplets is sufficiently small that droplet-droplet correlations can be 
ignored. The intermediate-field region, in which Eq. (|12D holds, we call the multi- droplet 
region. It corresponds to the picture of decay through continuously nucleating and growing 
droplets, ffist introduced by Kolmogorov, Johnson and Mehl, and Avrami over fifty years 
ago [|H52|. 

For stronger fields, a picture based on localized droplets no longer adequately describes 
the system, which in that case decays via long-wavelength, unstable modes reminiscent of 
spinodal decomposition [^. The crossover field separating the multi-droplet region from 
this strong-field region was termed "the mean-field spinodal" in (I). It is located at a field 

i^MFSP < ^K [|g 



54 



We expect that Re < Ro for all H < i^MFSP- A recent argument 
based on a transfer-matrix calculation indicates that the single-phase correlation length in 
the metastable phase is given by ^ms ~ T/(4J— 2/i) at low T ^1|. A rough estimate for the 
field at which droplet theory breaks down can be obtained by requiring that 2i?c = ^ms- This 
yields i?^MFSP <'2K and specifically ifMFSp(0.8Tc) ~0.33 with the critical-droplet diameter 
2Rc ~ 1.3, independent of system size. The dynamics in the strong- field region beyond 
-f^MFSP will not be discussed in detail here. 



B. Finite-size effects 



Next we study the effects produced by a finite system size L. For simplicity we consider 
hypercubic systems of volume L'^ with periodic boundary conditions. Effects of heteroge- 
neous nucleation at free boundaries are discussed in Ref. |^. In the large-L limit. 



L > i?o > ^c 



(13) 



we use an approximate argument by which the system is partitioned into {Lj RqY^I cells 
of volume Rq. Each cell decays from the metastable state to the equilibrium state in an 
independent Poisson process of rate RqI = tQ^. The volume fraction is then self-averaging, 
with (t(0s)) approximately equal to the infinite-system result of Eq. (|12D, and the relative 
standard deviation is 

'=- — wj) — "^^°/^)^^' (i'^) 

where p ~ 1 is the relative standard deviation of a single Poisson process. The regime 
characterized by r <C 1 was termed "the deterministic region" in (I). It is subdivided into 
the multi-droplet and strong-field subregions discussed above. 
For smaller L, so that 

/?o > ^ > i?c , (15) 

the random nucleation of a single critical droplet in a Poisson process of rate L'^I is the rate- 
determining step. This is followed by rapid growth, until this droplet occupies the entire 
system after an additional time much shorter than the average waiting time before a second 
droplet nucleates. Therefore, in this case the characteristic lifetime is 

{t{<P.))^{L'l{T,H)y' 

^ L~'^[A{T)]-^\H\-'^^+'hm^ (16) 

with r ^1 and only a weak dependence on the threshold (ps- This single-droplet region is 
part of "the stochastic region" observed in (I). 

The crossover between the deterministic and stochastic regimes is determined by the con- 
dition Loc-Rq- We identify the crossover field with the "dynamic spinodal field" introduced 
in (I), and in the limit H—>-0 we explicitly obtain from Eq. ( p,Ob| ) 

We emphasize that, although i^DSP vanishes as L— >-oo, the approach to zero is exceedingly 
slow, especially for d=3 and larger. Therefore, -^dsp t^^Y well be measurably different from 
zero for systems that are definitely macroscopic as far as their equilibrium properties are 
concerned. (As an illustration, increasing L from 100 to 10^° for d=3 decreases the leading 
term in ifosp only to approximately one-half of its original value!) Recent exact results 
53| , |5^ indicate that if the limit T — ^ is taken before L —* oo, then ifosp should remain 



nonzero. An estimate for the mean lifetime at -ffosp iii the infinite-L limit is obtained by 
setting 

{t at iJosp) oc to{T, Hdsp) " 



vo{T,Hy)sp) 
oc L/Hbsp ~ L(lnL)^/('^-^) , (18) 

where the nonuniversal, temperature dependent proportionality factors have been dropped. 
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Finally, we consider the limit 

i?o > i?c > ^ ■ (19) 

In this case the volume term can be neglected in Eq. (^, and the free-energy cost of a 
droplet occupying a volume fraction 0s = ^(i?)/L'^ is F((^s) =L'^~-^(/)^'^~^)/'^S(T), so that the 
first-passage time to a given (pg is independent of H and diverges exponentially with L'^~^. 
Since the dynamics in this region of extremely weak fields or extremely small systems is 
similar to that on the coexistence line, H=0, we call it the coexistence region. The crossover 
field between the coexistence and single-droplet regions, called "the thermodynamic spinodal 
point" (THSP) in (I), is determined for a given 0s by Qd{Rc/LY = (ps, which yields 



IV. NUMERICAL RESULTS 

The mean metastable lifetimes are illustrated in Fig. |^, which shows the average 
first-passage time to m=0.7 (0s ~ 0.13) for two-dimensional, square Ising systems with 
L=64,. . . , 720 at 0.8Tc. The simulations were performed using the Metropolis single-spin- 
flip MC algorithm with sequential updates, as discussed in Sec. |T|. The lifetimes, in units of 
MCS, are plotted on a logarithmic scale vs. 1/\H\. Each data point is averaged over n=100 
independent realizations, and the errors are calculated from the empirical standard deviation 
a in t{m=0.7) as o j ^Jn. For given |if | and sufficiently small L the lifetimes decrease with 
increasing L as (t(0s)) ocL~^, in agreement with Eq. ( |T6| ) for the single-droplet region and 
as shown in Fig. 3 of (I). For larger L they converge to an L-independent value, as expected 
from Eq. ([T^) for the multi-droplet region. This L~'^ behavior was previously remarked in 
Refs. [p!4| , p!5| , p0| , ^ , as was the presence of the L-independent regime in Refs. [p!4| , ^pl| . The 



size-independent data in the intermediate- and strong-field portions of Fig. |^ belong to the 
deterministic region, whereas the size- dependent data in the weak-field, long-time portion lie 



in the stochastic region. The droplet-theory results of Eqs. ( p.Oa| ), (0), and (^) show that 
in hoth the multi-droplet and the single-droplet regions the field dependence of the lifetime 
is determined by the nucleation rate I{T,H), so that the asymptotic slopes of ln(t(0s)) vs. 
1/\H\'^^^ in the two regimes are related by a factor (rf+l). Based on the four points for 
L=64 in the single-droplet region of Fig. |] (l/|-ff| >20) the effective slope was estimated by 
a least-squares fit to be 0.36(1). This slope is indicated by the dashed line above the data 
in Fig. 1^, whereas the dashed line below the data has 1/3 this slope. The agreement with 
droplet theory seems excellent. However, the measured effective slope is almost 30% larger 
than the theoretically predicted asymptotic slope, S(0.8Tc) = 0.278 840... , which is obtained 



by using the numerically exact §ij^ value E(0.8rc)/J = 2.648 81... in Eq. (^). 

The explanation for this apparent discrepancy lies in the power-law prefactors in 
Eqs. (|10a|) and (|16]), which yield the effective slopes 



A.»^ "'rr. =A|gr + A. (21) 



with \ = {b+c)/{d—l) and A = S(T) in the single-droplet region described by Eq. (|T6D, and 
X = {b+c+d)/{d'^—l) and A = E(T)/{d+l) in the mult i- droplet region described by Eqs. QlOaj) 
and (jg). 

In Fig. ^ we show two-point finite-difference estimates for Aeg, based on the average 
lifetimes in Fig. |^, for pairs of fields, \Hi\ > \H2\. These estimates are plotted vs. the field 
value that occurs in the finite-difference version of Eq. (0), \H\ = \HiH2\ In \Hi/H2\/\Hi — 
H2\. The error bars are calculated as 



O-A 



\ 



^1 \ , / 0-2 



(^i); W2). 



H1H2 



H1—H2 



where cxj and (tj) are the empirical standard deviation and mean for the lifetime at \Hi\, 
with i=l or 2. For clarity only two system sizes, L=128 and 720, were included in the figure. 
The other sizes give results similar to those shown. The solid straight lines correspond to 
Eq. (^ with h + c = 2 and the exact S(0.8Tc) = 0.278 840.... The data points for both 
system sizes cluster close to the lower of these two lines in the whole multi-droplet region. 
However, in the strong-field part of the region the deviations of the L=720 data from the 
line are considerably larger than the one-o"A error bars, leading to unacceptably small x^- 
probabilities in a weighted least-squares fit. These deviations may possibly indicate the 
presence of small but statistically significant corrections to the droplet-theory result. To 
obtain an acceptable fit for L=720 we therefore successively eliminated the data point with 
the largest value of |if | until the x^-probability Q stabilized at a reasonable value greater 
than 0.1. The resulting fit, which includes twelve points in the field interval 0.03<|iJ|<0.15, 
gives b + c= 1.95(8) and S(0.8Tc)=0. 282(8) with (5=0.41, consistent with the parameters 
used to draw the solid straight lines. 

Only the smaller system penetrates fully into the single-droplet region in the field range 
for which data could be obtained with a reasonable amount of computer time. As seen from 
the figure, the agreement between the L=128 MC data and the solid line is also good in the 
single-droplet region, even though the error bars are larger there than in the multi-droplet 
region. 

The dotted straight lines in the figure correspond to Eq. (^) with b + c = 3, the value 
expected from the discussion following Eq. (|^), and the exact 2(0. 8Tc). The agreement 
with the data is far inferior to the the solid lines. Since the field-theoretical prediction of 
6=1 has been well confirmed by several numerical methods [P5|- |^ , |45| -^ , we believe the 



disagreement between the predicted and observed values of b+c must be ascribed to c, the 
exponent giving the field dependence of the kinetic prefactor. (An alternative explanation, 
that the field dependence of the droplet growth rate, fo~ \H\, might be wrong, is probably 
ruled out since this only would affect the fit in the multi-droplet region.) We therefore 
conclude that our numerical results agree well with the predictions of droplet theory, but 
with c~l instead of c=2. We believe the unexpected value of c is an expression of the 
nonuniversality of the kinetic prefactor, and that it is a consequence of the sequential MC 
update scheme. This view is supported by the simulations with random updates, which we 
report at the end of this section. 

In the strong-field region, droplet theory is not expected to be applicable. For fields 
beyond ifMFSp(0.8Tc) ~0.33, Aefj indeed exhibits pronounced oscillations with H, the first 
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of which can be seen in the right-hand portion of Fig. (j^). Further study at these extreme 
fields is left for future work. 

The crossover between the multi-droplet and single-droplet regions is clearly seen in Fig. |^ 
as a jump in Aeg. However, it is difficult to determine the corresponding dynamic spinodal 
field Hdsp directly from this figure with sufficient accuracy to verify Eq. (p!7|). Instead, we 
consider the relative standard deviation r, which is shown in Fig. ^ on a logarithmic scale 
vs. 1/\H\. The error bars in the figure are estimated by standard error-propagation methods 
as 

r / 2(n-l) ^'/' 



1 + 



m- 



n 



In agreement with our predictions, r crosses over from the behavior described by Eq. (p!4D 
for large |iif| to r^l for smaller \H\. We take as our estimate for -ffnsp the field -ffi/2 for 
which r=0.5. This field is determined for each value of L from the crossing of a weighted 
least-squares fit to Inr in the linear region of Fig. ^ with the horizontal line r=0.5. The 
resulting estimates are shown in Fig. |^a as I/H1/2 vs. L on a logarithmic scale, with error 
bars estimated from those in Fig. ^ by standard error-propagation methods. The solid curve 
in Fig. ^a, which agrees well with the data, represents the crossover condition LocRq, where 
Rq is given by Eq. ( p.Ob|) with b + c = 2 in agreement with the fits in Fig. |^, and with the 
proportionality constant adjusted to approximately minimize the weighted sum of squares. 
For reference a dashed line with the asymptotic slope 3/H(0.8Tc), expected from Eq. (|T7|), 
is also included. It was obtained by using the same proportionality constant as in the solid- 
curve fit and setting the power-law prefactor equal to unity in Eq. ( |10b|) . The values of if 1/2 
for L=128 and 720 are also indicated by arrows in Fig. ^. The asymptotic L dependence 
of the lifetime at the DSP, given by Eq. (p!8D , is illustrated in Fig. ^d. For each L this 
lifetime was obtained by interpolation between the two closest field values bracketing H1/2 
for which simulations had been performed, and the uncertainty in the resulting estimate was 
obtained from those in Figs. ^ and |^ by standard error propagation. The numerical results 
are consistent with the analytical prediction. 

The crossover field between the single-droplet and coexistence regions, -^thsp of Eq. (pO|), 
is not easily observable in our calculation because of the smallness of the numerical factor 
(d-l)S(0.8Tc)/(2dA;BTms) =0.382 204... . Nevertheless, for L=64 and cutoff at m=0.9 (0s~ 
0.03) the resulting prediction is 1/Htusp ~28.5. As shown in Fig. ^ which corresponds to 
Fig. 1^ except that the cutoff is m=0.9, this estimate agrees well with the numerical data. 
A similar crossover was also observed in Ref. []K[. Except for the behavior for L=64 at low 
fields, and larger uncertainties in general. Fig. ^ is similar to Fig. 0. 

One can ask several questions regarding the effects of different local dynamics on the 
metastable lifetimes. The difference between the Metropolis dynamics used here and Glauber 
dynamics should be minor, at least for relatively weak fields [^, and we expect our main 



conclusions to remain valid for the latter as well. Perhaps more surprisingly, simulations 



with a modified Swendsen-Wang algorithm |T^j2yj5^ also seem to agree qualitatively with 
the droplet-theory predictions. 

The data presented above were obtained with a sequential-update algorithm. This is 
fast and well-suited to the architecture of the m-TIS2 computer, but it is clearly an unre- 
alistic depiction of any physical dynamics. To estimate the effects of the update scheme we 
also simulated the L=128 and 720 systems using the same Metropolis MC algorithm as in 
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the sequential-update study, but with updates at randomly chosen sites. This study was 
conducted entirely on heterogeneous workstation clusters. 

The most noticeable difference between the two update schemes is that it takes a con- 
siderably larger number of MCS to reach a given cutoff magnetization with random up- 
dates than with sequential ones. This effect and its nontrivial field dependence are il- 
lustrated in Fig. 1^, which shows the ratio between the random and sequential lifetimes, 
(tR(m=0.7))/(ts(m=0.7)), for L=128 and 720. Since the universal part of the field depen- 
dence must cancel in this ratio, its dominant field dependence is obtained from Eqs. (|10aj) 
and ( P^ as 

\^r(0s)/ I \H\ ''+1 in the multi-droplet region /nt^\ 

(ts(0s)) [ |iJ|-(cR-cs) in ^iie single-droplet region, 

with proportionality constants that depend only on T. Here cr and cs are the nonuniversal 
kinetic-prefactor exponents for random and sequential updates, respectively. This pure 
power-law behavior is only seen for L=720 in the weak-field portion of the multi-droplet 
region, indicating that additional nonuniversal effects distinguish between the behaviors 
of the two update schemes for stronger fields. Following the procedure described in the 
discussion of Fig. ||, we obtained a weighted least-squares fit to the L=720 data, including 
eight points in the field interval 0.032<|if |<0.09. This fit, which is shown as a solid line in 
the figure, yields (cr—cs) =1.24(3) with (5=0.25. In view of the very short fitting interval, 
this estimate should only be taken as an indication that the difference between the two 
values of c is probably on the order of unity. 

A somewhat puzzling feature in Fig. ^ is the fall-off of the ratio for L=128 in the single- 
droplet region. We believe this may be due to the proximity of the THSP for this relatively 
small system size. 

In Fig. p, which corresponds to and should be compared with Fig. |^, we show the effective 
slope Aeff from Eq. ( pT]) for the random-update case. The estimated -ffi/2 for both L=128 
and 720 (indicated by opposing arrows in the figure) agree to within the statistical error 
with the estimates for sequential updates. In the weak-field portion of the multi-droplet 
region the effective slopes follow the lower dotted line, which has the same meaning as in 
Fig. ^. Two major differences from the sequential-update case are apparent. Firstly, the 
field interval in which Aeff seems to follow a linear approach to the exact S(0.8Tc)/3 is 
narrower, as one would expect from the behavior of the lifetime ratios shown in Fig. |^. 
Secondly, the fact that the effective slopes cluster around the dotted line indicates that the 
nonuniversal prefactor exponent c is close to two, as expected for dynamics that can be 
described by a Fokker-Planck equation p4|-p^. Following the procedure described above. 



we obtained a weighted least-squares fit to the L=720 data, including nine points in the 
interval 0.032<|ilf |<0.11. The resulting parameter estimates were: S(0.8Tc)=0. 259(11) and 
c=2. 43(15) with Q=0.25. In view of the short fitting interval we consider that this estimate 
is consistent with the expected exact values. For |if |>0.11 the deviation of Agff from the 
dotted straight line takes the form of a smooth, slow oscillation whose amplitude is much 
larger than the statistical errors. 

In summary, our results indicate that the main difference in the nonuniversal part of 
the field dependence of the metastable lifetimes between the sequential and random update 
schemes lies in the kinetic-prefactor exponent c. For sequential updates we find c~l, whereas 
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for random updates c~2, in agreement with the theoretical expectation c=2. In addition to 
the difference in c for weak fields, the quantitative differences between the lifetimes in the 
two update schemes become progressively larger with increasing \H\. 

V. DISCUSSION 

We have demonstrated that the metastable lifetimes for impurity-free kinetic Ising mod- 
els with nearest-neighbor interactions and local dynamics, in weak to moderate unfavorable 
applied fields, exhibit system-size and field dependences that can be explained by a field- 
theoretical droplet model of homogeneous nucleation p3|-p6|. The most significant deter- 
mining factor for the lifetime is the free energy of a critical droplet of the stable phase, which 
in a wide field interval leads to an exponential dependence on the inverse applied field. This 
free energy is independent of the details of the local dynamics, and it can be accurately 
calculated by combining a Wulff construction with the exact zero-field equilibrium surface 
tension. 

Our numerical estimates for the lifetimes, based on Monte Carlo simulations in two 
dimensions, are in excellent agreement with the droplet model. This agreement includes both 
the aforementioned exponential dependence on the inverse field, and multiplicative power- 
law prefactors. There are two independent contributions to these prefactors: a universal 
exponent b, equal to unity for d=2, and a nonuniversal exponent c, which depends on the 
details of the dynamics. Our results indicate that with MC updates at randomly chosen 
sites, c is close to two, which is the exact result for dynamics described by a Fokker-Planck 
equation pi,|23|. With sequential updates we found c~l. It is tempting to conjecture that 



the exact value in this case is unity, but at present we have no theoretical arguments to 
substantiate this suggestion. For stronger fields we found additional, dynamics-dependent 
deviations from the droplet-theory predictions for Aeff. 

We have identified four different field regions, in which the decay proceeds through 
different excitations. In order of increasingly strong unfavorable field \H\, these were the 
"coexistence region", characterized by subcritical fluctuations on the scale of the system 
volume; the "single-droplet region", characterized by decay via a single critical droplet; 
the "multi-droplet region" , characterized by decay via a finite density of droplets; and the 
"strong-field region", in which the droplet picture is inappropriate. The crossover fields 
between these regions, 

[i^THSP-^"'] < [i/DSP~(lni:)"^] < [HuFSP ~ L"^] , (23) 

are accurately predicted by droplet theory. The different regions and crossover fields are 
illustrated in Fig. |I|. We believe the slow, logarithmic vanishing with system size of the 
dynamic spinodal field Hdsp, which separates the single-droplet and multi-droplet regions, 
is of particular significance. This convergence is so slow, especially in three and higher 
dimensions, that systems that are rightfully considered macroscopic in terms of their equi- 
librium properties may nevertheless possess a measurable single-droplet region in which the 
metastable state is exceedingly long-lived on the average. 
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FIGURES 

FIG. 1. Schematic sketch showing the field dependence of the average nietastable hfetime for a 
d-dimensional kinetic Ising model on a logarithmic scale. The dynamic spinodal point is indicated 
by the solid vertical line labelled DSP. To its left is the deterministic region, and to its right is the 
stochastic region. The deterministic region is divided into the strong-field (SF) and multi-droplet 
(MD) subregions by the mean-field spinodal, indicated by the dashed vertical line labelled MFSP. 
The stochastic region is divided into the single-droplet (SD) and coexistence (CE) subregions by 
the thermodynamic spinodal, indicated by the dashed vertical line labelled THSP. This figure is 
based on the analytical and numerical results presented in this paper, but it is not drawn to scale 
for a specific system size. In particular, the extent of the single-droplet region is greatly reduced, 
relative to the multi-droplet region. 

FIG. 2. Mean metastable lifetime for a two-dimensional kinetic Ising model at T=0.8Tc, es- 
timated as average first-passage times to m=0.7 {(ps ~ 0.13). The simulation was performed by 
the Metropolis MC algorithm with sequential updates. The lifetimes are given in units of Monte 
Carlo steps per spin (MCS) and are shown on a logarithmic scale vs. 1/\H\. The symbols indicate 
data for system sizes L=64 {+), 128 (Q), 256 (x), 400 (□), and 720 (O), ah obtained with a 
special-purpose m-TIS2 computer, and L=128 ({^) and 720 (filled diamonds), obtained on hetero- 
geneous workstation clusters. The solid curves are merely guides to the eye. The vertical arrow 
indicates the inverse field at which droplet theory is expected to break down, l/i?MFSP^3.0. The 
L-dependent data in the upper right-hand sector lie in the single-droplet region, the L-independent 
data in the central portion of the figure lie in the multi-droplet region, and the L-independent data 
in the lower left-hand sector belong to the strong-field region. The dashed line above the data has 
a slope of 0.36, estimated from the four points for L=64 with l/|i:/^| > 20, and the dashed line below 
the data has 1/3 this slope. The standard deviations in the average lifetimes range from on the 
order of the symbol size for the weakest fields shown, to two orders of magnitude smaller for the 



strongest fields, and error bars are therefore not shown. See details in Sec. IV. 



FIG. 3. Finite-difference estimates for the effective slope A^g from Eq. (pTD, shown vs. \H\. 
The parameters and symbols are the same as in Fig. ^. For clarity only data for L=128 (Q and 
0) and 720 (empty and filled diamonds) are included, and error bars are shown only where larger 
than the symbol size. The solid straight lines correspond to Eq. (p|) with b+c=2 and intersect the 
vertical axis at the exact H(0.8Tc) and H(0.8Tc)/3, respectively. The lower solid line gives a good 
fit to the data in the multi-droplet region for both system sizes. Only L=128 penetrates fully into 
the single-droplet region, where the agreement with the theoretical prediction (upper solid line) is 
also good. The dotted lines, which correspond to b+c = 3, do not fit the data as well. The crossover 
between the multi- and single-droplet regions is seen as a jump in Aeg . The vertical arrows with 
horizontal error-bar "feathers" indicate the positions of the estimator H1/2 for the crossover field 
-f^DSPi obtained from Fig. ^a. The left pair of arrows corresponds to L=720, and the right pair to 
L=128. The vertical arrow marked -ffMFSP marks the L-independent field at which droplet theory 



is expected to break down. See details in Sec. IV 
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FIG. 4. The relative standard deviation r for i(rn,=0.7), shown on a logarithmic scale vs. 
1/\H\. The parameters and symbols are the same as in Fig. |2|. The behavior of r crosses over 
from the approximate straight line described by Eq. (|l^ in the multi-droplet region to r ~ 1 in 
the single-droplet region. The inclined solid lines are weighted least-squares fits to the data in the 
region l/|i/|>10, r<0.6. Their average effective slope is 0.13(1), which is intermediate between 
H(0.8Tc) and H(0.8Tc)/3. The estimates -ffi/2 for -ffosp are found where these lines cross the 
horizontal dashed line at r=0.5. These estimates are also shown vs. L in Fig. |5|a and are indicated 



by arrows in Fig. y. See details in Sec. IV 



FIG. 5. Inverse field and mean lifetime at the DSP. (a): The estimates I/H1/2 for l/f^osp as 
obtained from Fig. |^, shown vs. L on a logarithmic scale. The solid curve corresponds to LocRq 
with b+c = 2 in Eq. (|lOb| ) and with the proportionality constant adjusted to minimize the weighted 
sum of squares. The dashed line indicates the asymptotic slope 3/H(0.8Tc). (b): The average 
lifetimes at I/H1/2. The solid straight line is a weighted least-squares fit. The figure supports the 
asymptotic behavior given in Eq. (|lH), {t at ^dsp) ^ LlnL. See details in Sec. |^ 

FIG. 6. Mean metastable lifetimes at r=0.8rc, estimated as average first-passage times to 
m=0.9 {(ps ~ 0.03). The parameters and symbols are otherwise the same as in Fig. |2|. Except 
for somewhat larger uncertainties and the behavior for L=64 at low fields, this figure is similar 
to Fig. ^. The guide-to-the-eye curves are dotted, except for L=64, which has been drawn solid 
to emphasize the behavior for this smallest system. It crosses over between the single-cluster and 
coexistence regions at a ^s-dependent -ffTHSP~28.5 (marked by the opposing arrows), as discussed 



in Sec. IV 



FIG. 7. The ratio {tfi)/{ts) between the mean lifetimes obtained with random and sequential 
MC updates for L=128 and 720. The symbols and parameters are the same as in Figs. ^ and|^. The 
two pairs of opposing arrows with horizontal error-bar "feathers" mark -ffi/2 for the two system 
sizes. The solid straight line in this log-log plot is a weighted least-squares fit to Eq. (p^) in the 



weak-field part of the multi-droplet region for L=720. See details in Sec. IV. 



FIG. 8. Finite-difference estimates for the effective slope Aefj from Eq. (p|) for random MC 
updates with L=128 and 720, shown vs. \H\. The figure is analogous to Fig. ^, and the symbols 
and parameters are the same as in that figure. The behavior is similar to the sequential-update 
case, except that the prefactor exponent 6-|-c~3, and the field interval in which Agfj seems to follow 



a linear approach to the exact S/3 is narrower than in Fig. |^. See details in Sec. |V 



